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Abstract. - We describe a formally exact method to calculate the optical conductivity of a 
one- dimensional system with fluctuating order. For classical phase fluctuations we explicitly 
determine the optical conductivity by solving two coupled Fokker-Planck equations numerically. 
Our results differ considerably from perturbation theory and in contrast to Gaussian order 
parameter fluctuations show a strong dependence on the correlation length. 



In a strictly one-dimensional system with only short-range interactions order parameter 
fluctuations prevail over long-range order. There is no phase transition at any finite temper- 
ature. Although strictly one-dimensional systems do not exist in nature a large number of 
organic and inorganic compounds can be considered to be electronically quasi-one-dimensional. 
Their conductivity in a pronounced direction can be orders of magnitude larger than in any 
other direction. An intensely studied class of such quasi-one-dimensional materials are Peierls 
systems. Below a critical mean field temperature mean field theory predicts long-range 
order but one-dimensional fluctuations prohibit the occurrence of this symmetry broken state. 
Only at temperatures much smaller than couplings between one-dimensional fluctuations 
become important and result in a flnite critical transition temperature Tc below which three- 
dimensional long-range order occurs as a charge density wave, a static lattice distortion, and 
a gap in the electronic spectrum Precursors of the ordered phase can be seen far above 
the Peierls transition temperature Tc (but below T^^) ^ where one-dimensional fluctuating 
order dominates. A strong suppression of the density of states, i.e. a pseudogap, is quite 
similar to the pseudogap found in the cuprate superconductors and can also be observed in a 
number of correlation functions. 

The instability of a Peierls system is due to the coupling of its conduction electrons to 
quasi-one-dimensional acoustic phonons. At not too small temperatures it is usually assumed 
that the phonons can be approximated by a static order parameter fleld A(a;) which for 
an incommensurate lattice filling turns out to be complex. Above the phase transition the 
order parameter field can be considered to resemble disorder with short-range order such that 
its expectation value {A(x)) has to be zero. The second moment decays exponentially, i.e. 
{A{x)A*{x')) = A^e^'"^^^ While the density of states of an incommensurate system is 



2 



essentially unaffected by white noise 0, but D = A^^ — const) a pseudogap is 

established for sufficiently large correlation lengths ^ [3ini|71|l|HlEl[Tni[niiniiniIIlIISlIini- 

Numerical simulations show that for Gaussian statistics the density of states vanishes as 
p{0) oc (with /i w 0.64) at the Fermi level ^HI- For classical phase fluctuations, however, 
it can be shown analytically that p{0) oc exp(— 4As^) jT^- Hence, the finiteness of the 
correlation length and even the statistics of higher moments turn out to be very important. It 
is therefore obvious that a non-perturbative method which is not restricted to Gaussian colored 
noise should also be needed for an accurate calculation of one- or two-particle correlation 
functions. 

Experimentally measurable transport properties can be related to two-particle correlation 
functions. Of particular interest is the optical conductivity which is proportional to the optical 
absorption of a solid. For a one-dimensional system with Gaussian white noise the optical 
conductivity is given by the generally accepted Mott-Berezinskii law [17iil8u 3u4iil9 ) 



Here ctq is the Drude conductivity, lo is the frequency and r is the elastic scattering time. Of 
course it would be interesting to see how the Mott-Berezinskii law is modified by rendering the 
correlation length finite, opening also the possibility of studying non-Gaussian fluctuations. 
The optical behavior of a system with fluctuating order would also be interesting within the 
light of optical investigations of the pseudogap phase of quasi-one-dimensional conductors like 
(TaSe4)2l or K0.3M0O3 above their charge density wave transition |20inil21ll22| . The origin of 
the pseudogap observed in the optical conductivity of these systems can be considered to be 
complementary to the origin of a gap in a one-dimensional Mott insulator where the optical 
gap is dynamically generated by electron-electron interactions |23[I24) . 

In this paper we will show how Halperin's method to calculate correlation functions for a 
particle in a white noise potential j25l4| can be generalized to calculate the optical conductivity 
of the fluctuating gap model which encapsulates the essential features of the pseudogap phase 
of quasi-one-dimensional Peierls systems and is also of interest for other systems [5_ 26 6 7 4I14| . 
The fluctuating gap model describes electrons with a linearized energy dispersion interacting 
with a static space-dependent order parameter field A(x) scattering electrons from one Fermi 
point to the other. The non-perturbative calculation of the spectral function for this model 
with disorder determined by a generalized Ginzburg-Landau equation was addressed in a 
previous publication jl6| . 

As the order parameter field A(a;) is complex for an incommensurate lattice filling, fluc- 
tuations of A(x) can be subdivided into fluctuations of its amplitude A{x) and its phase 
'&{x), such that A{x) = A{x)e^^'^^\ Using natural units which make Planck's constant h, 
Boltzmann's constant ks, and the Fermi velocity vp equal to unity the Hamiltonian of the 
fluctuating gap model is given by 



Since H is invariant with respect to the charge conjugation operation H —>■ cTiH*ai where ui 
is a Pauli matrix we can apply boundary conditions which imply that all eigenstates "^nix) of 
H are simultaneous even eigenstates of the charge conjugation operation xpn{x) (Jixp*^{x) 

such that Ipnix) = {ll^nix),tpnix))'^ ■ 

In the following we will be interested in the optical conductivity cri (lu) at flnite temperature 
T which can be expressed by the Kubo formula 



(1) 




(2) 
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ai{uj) = vrWp 



de 



F(e,e + cj) . 



(3) 
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Here tOp = Airne^ /ra* is the square of the plasma frequency (with — e and to* the charge and 
effective mass of an electron and n the density of electrons), /(e) = [exp(e/T) + 1]^^ is the 
Fermi function and F is defined as 

^(^'^') ^ Z \Jnn.?^{e - e„)J(e' - e™)\ , (4) 

\n,m / 

where 

Jnra = j dx Xpl^i^x) (J m{x) (5) 



is a current matrix element. The process of averaging in Eq. Q restores translational sym- 
metry. It is therefore possible to rewrite Eq. Q as 

F(e,e') = (E / dxV:(x)V™(a:)V;„(0)Vn(0)<5(e-e„)<5(e'-e^)\ . (6) 

Let US now decompose the wave function i)n{x) according to 

{x) = lAe^'r. (-)+*k.„ (-)+'?(-)]/2 . (7) 

It should be noted that t^^^ [x) measures the phase of 'ipn{x) relative to the phase of the order 
parameter field i?(x). From Hxpn{x) = e„i/'„(a;) it follows immediately that the phase variables 
and satisfy the following equations of motion 14 , 

d^ipeix) = 2e~ d^-d{x) +2 A{x) COS [ip^{x)] , (8) 
d^x) = A{x)sm[ip,{x)] . (9) 

To construct all eigenstates for a given realization of the disorder potential we follow Ref. ^ 
and integrate the equations of motion (|8I9|I from the boundaries at x = tL/2 to x = 
where we try to match these solutions. If Lpe,±{x) and (e^±{x) are the solutions to Eqs. (|8I9(I 
when advancing the solution from x = ^L/2 towards x = 0, the matching condition which 
determines the discrete set of energies e„ corresponding to the eigenstates ipn{x) is 

^e,+ (0) = ^e,-(0) . (10) 

Instead of working with normahzed wave functions which satisfy / dx ilil^{x)il'm{x) — Sn^m let 
us from now on use the normalization condition !■!/)„ (0)| = 1 which amounts to the replacement 



^"^""'^ i2jdx\Mx')\^y/^ ■ ^ ' 



Making use of the identity |H[T^ 



^(^,,(0) - ^e,-(o)) = E ,rJT::'.. ■ (12) 

" 2j_^/2dx'|V'„(a;')r 

which can be shown by using the theorem 5{f{x)) — \ f\x) \ '^(^ ^ '^here the 
the zeros of /, we can now rewrite F as 

F(e,e')-2Re^W+(0)<5(^,.+(0)-^,,_(0))%,,,+ (0)-^,.,_(0))^ . (13) 
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The quantities W±{x) are defined as 



W±{x) 



and satisfy tlie equation of motion 

d^W{x) ^-i[e' -e + A{x)e^''^'' - A{x)e"<''] W{x) + 1 



(14) 



(15) 



At this point let us specify the process of averaging. We assume that the probabil- 
ity distribution of {A(a;)} can be expressed via a generalized Ginzburg-Landau theory, i.e. 
P{A} = ie-^{^>/^, where 



F{A}/T = 



L/2 



dx 



L/2 



id.Aix)r + ^ id.,ix)r + viAix)) 



(16) 



M{A), n{A), and V{A) can in principle be evaluated microscopically but can also be taken to 
be phenomenological given functions. At low temperatures, amplitude fluctuations get frozen 
out, i.e. A{x) « As and Eq. Hl()|l reduces to a form describing classical phase fluctuations, 



L/2 

Fm = ^ I dx {d,^{x)f 

L/2 



(17) 



Here Us = Tn measures the stiffness of long- wavelength fluctuations and is approximately 
independent of temperature. A two-dimensional analogue to Eq. (|17|l has been used by Emery 
and Kivelson '57' to explain the pseudogap phase in the underdoped cuprates. 

The essential property of Eq. ((TBI) is that the stochastic process corresponding to Eq. l(TBjl 
has the Markov property and hence can be described by the Langevin equations 



d^A{x) = a{A)+b{A)r]A{x) , 
d.,^{x) = n-^/^{A)7ji,{x) , 



(18) 
(19) 



where r]i{x), i = A, i9 are independent Gaussian white noise stochastic processes with {rii{x)rii{x')) 
5{x — x'). The functions a(A) and h{A) may be obtained from Eq. (|16|l using the transfer 
matrix method |28[ll5l[TB| . Decoupling the average in Eq. (|13|l by introducing the auxiliary 
variable A we obtain 



F(e,e') = 2Re / dA / d^p dtp' 



,Pi (.^,y',A;0)Po(-v7,-y',A;0) 
P(A) 



(20) 



where P(A) = {5{A - A{x))) and 

P„(A,(^,^';a;) = ((W^+(a;))" (5(A - A(x)) ^(^ - ^,,+ (a;)) - ^e',+ (x))) , n = 0,1(21) 

are joint probabilities. Differentiating Eq. (|21|l with respect to x we obtain the following 
coupled three-dimensional Fokker-Planck equations 



2n 
2n 



- 2d^{e + Acosip) - 2d^'{e' 4- A cos (^') 

- 2(e -I- A cos (p)d^ - 2{e' + A cos (p')d^i 



Po , (22) 



z(e' 



Ae 



dAa 



Pi + Po .(23) 
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Fig. 1 - Optical conductivity calculated for classical phase fluctuations and T = 0. The correlation 
length is given by AsS, = 0.5, f , 2, 4, 8, and oo. Dashed line: extrapolation between cri(O) — and our 
numerical data for As(, — 0.5. 



The solutions to these equations determine the exact optical conductivity for our problem 
and are similar to Halperin's equations to determine the optical conductivity of a particle in a 
white noise potential 12514) . Unfortunately solving for the stationary solutions to Eqs. (|22I23|I 
is not an easy task. Both equations are parabolic partial differential equations. The absence 
of a diffusion term (d^p — d^'Y' makes their numerical solutions notoriously unstable. However, 
adding a small but finite diffusion term K{d^ — d^p'Y turns the equations into elliptic partial 
differential equations which are numerically accessible. To further simplify our problem we 
will restrict ourselves to the physically especially interesting case of classical phase fluctuations 
where 2n = ^ and A drops out of all equations. This reduces the dimensionality of the partial 
differential equations to two. Since the coefficients of both equations involve only exponentials 
of (/? and (/?' and since both equations are subject to periodic boundary conditions the Fourier 
transform method seems to be the method of choice to solve for Pq and P\ . We find that as 
wc increase the number of Fourier components, the parameter k can be chosen smaller and 
smaller to obtain a stable solution. For sufficiently small k our solutions become independent 
of K. 

In Fig. ^ we show the optical conductivity calculated for various correlation lengths and 
zero temperature. Unfortunately we have not been able to calculate oiiiS) for to < O.lAs. The 
difficulty to calculate u\{lS) for small uj can be understood as follows: As e' approaches e the 
phase i^iix) will approach V3g/(x) such that /o(</'j ~* The singularity in this 

probability distribution obviously causes numerical difficulties. As can be seen in Fig. ^ the 
restriction to frequencies w > 0.1 does only imply that we cannot use our numerical methods 
to determine the quantitative corrections to the Mott-Berezinskii law. All curves aim towards 
(Ti(0) = which is consistent with the fact that the dc conductivity of a one-dimensional 
disordered system should be zero. In the limit ^ — > cxd the optical conductivity approaches 
the well-known result of a Peierls insulator j26j . 

involving a square root singularity at w « 2A3. The asymptotic l/w'^-decay applies to all 
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Fig. 2 - Optical conductivity calculated for classical phase fluctuations, Aa^ = 0.5, f,2,4, 8, oo and 
temperature given by T = f/vr^. Dashed lines: extrapolation between a"i(0) = and our numerical 
data. 

finite correlation lengths ^ and shows that for large oj perturbation theory is accurate. In 
the pseudogap regime, however, ai{uj) shows a strong dependence on the correlation length 
^. This is in contrast to the case of Gaussian order parameter fluctuations (for which there 
exists an "almost exact" solution that has also been applied to a two-dimensional model of the 
pseudogap state, see Ref. where the dependence on ^ is rather weak and no singularity 
arises for — > cxd |^. Finally, as a nontrivial check for the optical conductivity we have the 
conductivity sum rule |29| 



which is well satisfied by our numerics. 

In the pseudogap phase of incommensurate Peierls systems the correlation length ^ is 
approximately given by ^ = I/ttT [T1|13| . We therefore show in Fig.|21the optical conductivity 
for different values of the correlation length and corresponding temperatures. Due to thermal 
excitations, weight of the peak at w w 2As is moved to lower frequencies. In this case the 
behavior of cti(w) for small uj is not so obvious. However, we know that all eigenstates of 
our strictly-one-dimensional system are localized and therefore the dc conductivity should 
vanish for any finite temperature 0]. This enables us to find an extrapolation of cfi{uj) in the 
frequency regime < w < 0.1. We would like to note that our results differ from perturbative 
results of the optical conductivity which were calculated within the same model to first order 
in perturbation theory without considering vertex corrections ^21- The perturbative results 
give a finite dc conductivity which for this model must be an artifact of perturbation theory 
but nevertheless explain the qualitative behavior of the dc conductivity of Peierls systems 
above their charge density wave transition. 

In real Peierls systems dynamical effects become important below a crossover energy scale 
a;* <IC Afi. If pinning is only weak, a sliding charge density wave can form and lead to an 
enhanced optical conductivity at low frequencies. This collective dynamic behavior cannot 
be described by our model. Above the energy scale w*, however, our model is applicable and 
describes the peak structure at w « 2As. 




(25) 
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In summary, we have presented an exact method based on the phase formahsm to calcu- 
late the optical conductivity of a one-dimensional system with fluctuating order. The optical 
conductivity can be expressed through the solutions to two coupled partial differential equa- 
tions which for classical phase fluctuations we have solved numerically. Our results correct 
previous perturbative results and give an explanation of the precursors of a single particle gap 
in the optical conductivity. To explain the experimentally observed collective features at low 
frequencies the dynamics of the order parameter field should be included. It would also be 
interesting to study corrections to the Mott-Berezinskii law for finite correlation lengths and 
different statistics of higher moments of the order parameter field. The governing equations 
for such a calculation are our Fokker-Planck equations. 

* * * 

I would like to thank Peter Kopietz and Hartmut Monien for helpful discussions. 
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